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1.0 SUMMARY 


The radar baekscatcer from the sea depends upon both the magnitude 
and direction of the wind vector. This dependence can be explained by 
both theoretical and expiries! node Is. These results are extresely 
important since in the near future NASA proposes to launch the SeaSat-A 
oceanographic satellite. One of the goals of this satellite is to demon- 
strate the feasibility of global wind speed and direction measurements 
over the ocean. Such wind vector information will provide meteorologists 
with essential input into wave and weather forecast models. 

Several regression models are proposed to explain the wind direction 
dependence of the 1975 JGMSWAP (Joint North Sea Vfave Project) scatter©- 
meter data, in which continuous directional observations were made by 
flying in a circle. The models consider the radar backscatter as a 
harmonic function of wind direction. The constant term accounts for the 
major effect of wind speed and the sinusoidal terms for the effects of 
direction. The fundamental accounts for the difference in upwind and 
downwind returns, while the second harmonic explains the upwind-crosswind 
difference. Possible higher order effects would appear as coefficients 
of the higher harmonics. The first model considers only terms involving 
the fundamental and second harmonic. The second and third models are 
more general, since third and fourth harmonics are included, respectively 

The results, for angles of incidence of k0° and 65°, indicate that 
third order effects are negligible. Also, for all cases, the first model 
which used only cosine terns up through the second harmonic, appears to 
explain the angular variation adequately. The coefficients of this 
model are related to the wind speed by a power-law expression. The rela- 
tionship is extremely valuable, since it simplifies the process needed 
to determine a wind direction from two orthogonal scatterometer measure- 
ments. Thus, the results are of value to future satellite observations 
of the ocean. 


An inversion tecEcncque is also prcjrosed which will eliainatte the 
effects of wind speed and direction. This approach assies that two 
orthogonai scattering measurements are available. Also, if swell or 
the large-scale waves do not line op with the wind raves, it is asserted 
that the necessary information is available to correct for these effects. 
The technique performs adequately under these constraints, and also, 
when possible higher order wave effects are present. 

A comparison between Two different types of sea scatter theories, one 
type represented by the work of Wright and the other by that of Chan and 
Fung, was made with recent scattercmeter measurements. Wright's theore- 
tical model differs from that of Chan and Fung in two major aspects: 

(I) Wright uses Phi 51 ips 1 model for the sea spectrum, while Chan and Fung 
use Hitsuyasu and Honda's spectrum as proposed by Pierson; (2) Wright uses 
the modified slick-sea slope distribution of Cox and Hunk defined with 
respect to the horizon, while Chan and Fung use Cox and Hunk's slick-sea 
slope distribution defined with respect to the plane perpendicular to 
the look direction. 

Satisfactory agreement between theory and experimental data were 
obtained when Chan and Fung's model was used to explain the wind and 
azimuthal variation of the radar backscacter. Wright's theory does 
not include anistropic characteristics, and hence, can be only used 
to show wind speed and incidence angle dependence. Since Phillip's 
sea spectrum is independent of wind speed, Wright's modei does not 
show significant wind dependence, and thus, leads to poorer agree- 
ments with experimental data. 
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2.0 AM EMPIRICAL MODEL FOR OCEAN RADAR BACKSCATTER AMD iTS APPLICATION 
IN INVERSION ROUTINE TO ELIMINATE WIND SPEED AMD DIRECTION EFFECTS 


A. Modelling the Wind Speed Dependence of Sea Return 

The ocean surface can be viewed sirapl istically as consisting of 
waves of two scales: (i) gravity waves which have wavelengths ouch 

longer than those of the radar, and (ii) capillary waves which have 
wavelengths on the order of or smaller than the radar wavelength. At 
microwave frequencies and for incidence angles in the range 80“ > flj 1 25°, 
the radar return, especially the vertically polarized case, can be ade- 
quately explained in terns of the first order Bragg scattering theory. 

The selected ocean wavelength predicted by this theory is related to 
the radar wavelength by 

Aslna*~ (2-1) 

where A = ocean wavelength 
A = radar wave length 

Thus, a strong dependence of the radar return upon the sea spectrum in 
the capillary region is expected. Pierson (1975), Mitsuyasu and Honda 
097*0, and others have shown that the sea spectrum in the capillary 
region grows with increasing wind speed. For this reason, a strong 
correlation between radar sea return and the wind speed over the sea 
surface is anticipated. 

Work by many experimenters, most recently Young and Moore (1976), 
and Jones v Schroeder, and Mitchell (1977) suggest that the wind speed 
can be determined from the radar backscatter by using an empirical rela- 
tion of the form: 

o° = AU y (2-2a) 

or 

aJ B = A dB + Y(10 log U) (2-2b) 
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where A = levelling constant 
U = wind speed 

Y e* wind speed power coefficient. 

The effect of wind direction on sea returns has not been as 
extensively explored as the magnitude dependence. The majority of 
research in this area has been conducted under the NASA AAFE (Advanced 
Aircraft Flight Experiment) program. Attempts at modelling this effect, 
however , have never been truly satisfactory, since a relationship between 
the coefficients of the models and wind speed has not been shewn pre- 
viously (Afararci 0975))* 

Several regression models were applied to the 1975 J0H5VAP (Joint 
North Sea Wave Projected rcle-f I ight data in an attempt to determine 
the wind direction and speed dependence of the radar returns. The pro- 
posed models assume the hackscatter to be a harmonic function of wind 
direction. The constant tern accounts for the major effects of the wind 
speed magnitude and the sinusoidal terms account for the effects of wind 
direction. The fundamental accounts for the difference in the upwind 
and downwind returns, while the second harmonic explains the upwind-cross- 
wind difference. Possible higher order effects which nay occur are 
accounted for in some of the models by the higher order harmonics. To 
include these possibi 1 ities the following regression models were 
proposed. ^ 

Model 1 : o°(u,0,«) A^tU.e) cos' fnte-.J.fl (2-3a) 

3 

Model 2 : a°(u,0,$) *2L A n (u,9) cos Jn($-$ u )| (2-3b) 

k 

Model 3 : o°(U,8,$) =2 A (U # 0) cos fn ($-$)! (2-3c) 

n=0 L J 

where A n (U,6) = coefficient to be determined through regression 
$ - measurement direction from North, in degrees 
<fr u ~ upwind direction, in degrees 
U - wind speed, m/s 

8 » incidence angle (relative to vertical) 

*( Jones, Schroeder, Mitchell, 1978) 

-U- 


2 

The multiple correlation coefficient, R , was chosen to provide 

2 

a numerical evaluation of the regression analysis. R is a measure of 
the “proportion of total variation about the mean explained by the regres- 
sioni" a perfect fit has a value of one (Draper and Smith, 1966). 

By fitting Model 3 to the measured data, values tor the A n *s are 
obtained as shown in Table I; These A n 's are for a given incidence angle, 
6 , 'and a given flight (or wind speed). Reasonable improvements in the 
value due to use of third and/or fourth harmonics would imply that 
higher order effects are significant. As Table l shows, this is not the 
case. Even flight ?6, which has unequal crosswind values, appears to 
be adequately explained by Model 1, as shewn in Figure I. 

Since the second-harmonic model adequately explains the wind direc- 
tion effects on backscatter, an attempt was made to determine a rela- 
tionship between the three regression coefficients of this model and wind 
speed. The scattering coefficient is known to have a power-law dependence 
on the wind speed. A power-law relation of the following form is assumed 
for the A n 's, 


. Y_( 0 ) 

A n (u,e) = p n (e) u n 


(2-4) 


where p^fe) - levelling coefficient 

y (9) =* wind speed paver coefficient 
This expression is fitted to the A n 's shown in Table 1 at a given 
incidence angle. Figure 2 is an illustration of this case when 0 - 40°. 

It is clear from Figure 2 that the assumed expression for A n does provide 
satisfactory fits, with the possible exception of Aj (VV) . Table 2 shows the 
coefficients in (4) obtained through regression analysis for 8 = 40°(Figure 2). 
This relationship is also shown using both AAFE and JONSWAP results in 
Appendix B. 

B. VJiHDVEC - A Simple Inversion Technique to Eliminate Wind Speed and 

Direction Effects from SeaSat Orthogonal Radar Backscatter Measurements 
In the near future (1978), NASA intends to launch the SeaSat-A 
oceanographic satellite, carrying as one of its instruments the SeaSat 
Active Scatterometer System (SASS). The SASS will use four antennas 
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Aspect Angle (Degrees) 


Figure la: Rcf|t*oc>r.lon and observed values of plotted 

versus the azimuth angle, U 8.9 m/s. 







HH Polarization 


Date (Flight 


8/29/75 

9/02/75 

9/03/75 

9/09/75 

9/09/75 

9/10/75 


Wind 

Speed 


JONSWAP Regression Results (o * 40°) 
(Individual Flights) 


1.91 

2.01 

4.74 

15.41 

11.13 

5.96 


0.58 

1.05 

2.02 

5.07 

3.71 

2.02 


0.85 

0.83 

2.20 

7.80 

4.93 

2.68 


0.63 

kkk 

0.54 


0.91 

0.98 

0.95 

0.97 

0.97 

0.95 


(b). VV Polarization 


Date 


* j 


9/08/75 

9/09/75 

9/09/75 

9/10/75 


Flight 

Wind 

Speed 

■Ell 

4.5 

HI 

5.5 

16 

8.8 

17 

12.8 

18 

11.3 

19 

7.5 


2.95 

4.00 

9.32 

27.31 

21.63 

11.60 


A u 

(xio 3 ) 

A 2_, 
(xIO 3 ) 

kkk 

1.62 

1.07 

1.90 

1.95 

5.43 

3.73 

15.54 

4.12 

12.86 

0.70 ; 

6.55 


k*k 

1.50 


kkk 

*** 

irk* 

•kirk 

-0.69 

kkk 


0.92 

0.94 

0.92 

0.97 

0.98 

0.97 


kkk 

kkk 

kitk 

0.98 

0.98 

0.98 


These terms of the regression were not statistically significant. 























JONSMAP Regression Results (e ■ 65°) 


(a). HH Polarization 


Date 

Flight 

1 9/02/75 
9/0S/75 
9/09/75 
9/09/75 

14 

16 

17 

18 



5.5 

8.9 

12.3 

10.5 


0.17 

0.29 

1.26 

0.97 


0.068 

0.18 

0.83 

0.69 



0.032 

0.16 

0.55 

0.47 


0.84 

0.93 

0.96 

0.92 


0.86 

0.94 

*** 


(b). VV Polarization 


Date Flight 


9/02/75 

9/08/75 

9/09/75 

9/09/75 

9/10/75 



Wind 

Speed 

A 0 _ 3 
(xlQ 3 ) 

A 1 

<xto" 3 } 

A 2 

(xlO 3 ) 

5.5 

0.85 

0.098 

0.44 

8.9 

2.83 

0.89 

2.22 

12.3 

7.14 

1.70 

4.33 

10.5 

5.36 

2.75 

3.82 

7.5 

2.86 

t 

0.40 

2.22 



A 

(i 


0.28 

ft** 

0.30 


0.95 

0.93 

0.92 

0.83 

0.96 


ft** 

0.95 

0.94 

ftftft 

0.98 


*** 

ft** 

ft** 

0.99 


Flight 13 contained no 65° data. Also, the horizontal polarization data for flight 19 was deleted due 
to small number observations. 


These terms of the regression were not statistically significant. 
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Table 2: POWER LAW RELATIONSHIP OF REGRESSION COEFFICIENTS 

OF 2nd HARMONIC MODEL AT 0 = 40°. 


Coefficient 

Horizontal Polari; 

91 

Vertical i 

3 olariz< 

ition 

P N - 
xlO" 3 ) 

•V 

R 2 

|^1 

v n 

R 2 

A 

0 

7.53 

2.05 

0.94 

11.75 

2.13 

0.96 

A i 

3.46 

1.94 

0.96 

2.68 

1.95 

0.73 

A 2 

2.55 

2.16 

0.93 

5.02 

2.26 

0.95 























squinted 45° off the satellite subtrack to achieve "a star-like illumina- 
tion on the earth's surface," Grantham, et al.(1975). An operational 
frequency of 14.6 GHz was chosen to provide maximum wind sensitivity 
and minimal atmospheric absorption. Both vertical and horizontal like - 
polarization measurements will be made over an incidence angle range of 
[25°, 55°]. Doppler filtering will be used to divide this 1000 km. 
swath into either 25 km. or 50 km. cells as specified by the SeaSat 
User Working Group (UWG). 

The User Working Group has also specified that the SASS should be 
capable of measuring wind velocity within an accuracy of ±2m/s (U 20m/s) 
or ±10% (U 20m/s) for wind speed, and an accuracy for wind direction of 
±20°. This creates the need for an inversion technique, which wi 1 l use 
the orthogonal measurements, to infer these quantities. This problem can 
be solved in two steps: (1) wind speed, and (2) wind direction. 

1 . Wi nd Speed 

The wind speed can be estimated within the user's requirements by 
using the average for the same point from orthogonal viewing angles. 

This average is approximately equal to the radar cross section when the 
aspect angle is 45°, i.e. 


o°(<|>= 45°) 


AV 


oV*i> * °V*i + 9QO) 

2 


(2-5) 


From Model 1, it can be shown that (5) becomes exact for $ j = -45 and 
it provides the worst estimate when $ j = 135°. This is true for any 
incidence angle for which Model 1 is valid. It should also be noted 
that (5) provides poorer estimates when the upwind-downwind ratio becomes 
larger. 

As an illustration, the relationship between o°(iji = 45°) and wind 
speed is shown in Figure 3 for e = 40°. Again a power-law dependence 
of the form 

Yt (6) 

o°fo - 45°) - ^(eJu 1 (2-6) 
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Wind Speed (m/s) 


Figure 3: Wind speed power-law dependence of the 

radar cross section at $ - 45°. 
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provides satisfactory fits to the data. Consequently, the wind speed 
can be estimated knowing the average of the two scatterometer measure-* 
ments, i.e. 



where a = i/a,: a(HH) = 85.1, a(VV) * 60.Q 

y - 1/y,: y(HH) = 0.47, v(W) = 0.45 

As a test, wind speed errors were obtained from scattering coefficients 

calculated using the wind direction models. As seen from the results in 

Figure 4, the wind speed error can be expressed by 

U ^ = C + C. cos(^ + 45°). 
err o I i 

A wind speed correction of this form would offer considerable improvement 
if the data set had no fading or other random fluctuations. However, the 
SASS measurements will include fading and instrument bias. 

At present, no wind speed correction is employed. When this technique 
is applied to orthogonal pairs of JONSWAP measurements the results are still 

within the UWG's guideline of ±2 m/s (u <_ 20 m/s) or ±10? (u j> 20 m/s). 

As an illustration, the wind speed errors for flight 17 are shown in Figure 
5. Appendix C shows the wind speed error for the remaining JONSWAP flights. 

2. Wind Direction 

In the previous section, it was shown that the wind direction effects 
on o° can be modelled by expression (3a)- The regression coefficients of (3a) 
have a wind speed dependence given by Equation (4). In Section A, a method 

A 

was discussed for determining an estimate of the wind speed, U. Using 
this estimate of wind speed, the scattering measurement can be modelled by 


o 


o 


- A 2 cos2$ w + Ajcost^ + A q 


A A 

where A^ are estimates obtained from (4) using the U from (7). 
Consequently, the wind direction, <}> w , can be found from 


* w " cos 


-1 


-A, Wg - 8^(fl o - A, 

4a„ 


o°) 


(2-8) 
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(2-9) 


WIND SPEED ERROR Cm/s) 


Figure k&i Predicted wind speed errors vs U for HH-pofar izat ion 
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Figure hbt Predicted wind speed errors vs U for W-polarizatiop 
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Figure 5a: Wind speed errors vs U for flight 17, HH, 


i»0 * 


•SO 


0*0 * 


■*50 


* 1*9 


-1.5 ♦ 


*2*0 


•2*5 


-3.0 


-3*6 


Flight 17 0 
Incidonco Anglo t 40 
Polarisations HIS 


l 1 


1 . 


1 11 


i 1 

i i 


i i 


ii i • 

i t 

a 

t s 


+ 

V 


♦ 

♦ 


•*«*e ♦ 

•26c * 20c 60 o 100 i«f0 ISC 220 £60 SQO 3<»0 3S0 

0,8 «tl), 50 , t20 J60 200 2S0 . £40 320 360 

ASPECT ANGLE (Degrees) 


* 


J . 


VINO SPEED ERROR (ra/s) 


Figure 5b: Wind speed errors vs U for flight 1 7 , VV 
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Equation (9) holds for all incidence angles for which Node) 1 is valid. 

This inversion yields two estimated wind directions for each scattero- 
meter measurement, one in the first quadrant and one in the second. These 
are mirrored to obtain estimates in the third and fourth quadrants, 
respectively. Ho information is known from the measurements themselves 
to determine the wind quadrant, but in an operational satellite system, 
such as SeaSat, this can be obtained from V-IR or conventional meteoro- 
logical data. 

Examples of the wind direction errors obtained when applying this 
technique to JONSWAP flight 17 are shown in Figure 6. where A and B 
are Vhe fore and aft beam results, respectively. Considerable improvement 
can be obtained by averaging the results of both fore and aft measurements 
after they have been corrected to North. Results for flight 17 using this 
correction are shown in Figure 7» The wind direction errors produced by 
this method for vertical polarization are considerably lower than those 
for hbrizontal polarization, as seen by Table 3* One explanation for this 
observation is that the upwind-downwind scattering coefficient ratio 
increases with the wind speed much faster for the horizontally polarized 
case. Also, the scattering coefficient for the horizontal polarization 
can be much smaller than that for the vertical polarization and it drops 
off rapidly at large incidence angles. These later properties make it 
much more difficult to obtain reliable measurements. 


C. Conclusion 

The wind direction properties of radar backscatter were modelled 
by cosine fourier series through the fourth harmonic in wind direction 
from upwind. A comparison with JONSWAP data indicated that the effects 
of the third and fourth harmonic are neglible. 

A simple inversion technique, which uses two orthogonal scattering 
measurements, is also shown which eliminates the effect of wind speed 
and direction. The procedure used is as follows: 

(!) Estimate wind speed from average o° from 2 antennas 

(2) Estimate coefficients of directional equation by using o° 
estimate from (!) and previously determined wind-speed 
dependence of coefficients. 
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Figure 6a: Wind direction errors vs 0 using only 1 beam, HH 
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Figure 6b: Wind direction errors vs <J> using only 1 beam, VV. 
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Figure 7a : Wind direction errors vs $ using corrected averaged 

results of fore and aft beams, HH. 
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Figure 7b: Wind direction errors vs <}> using corrected averaged 

results of fore and aft beams, VV. 
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Table 3: WIND DIRECTION ESTIMATE ERRORS (IN DEGREES) 


FLIGHT 

HORIZONTAL POLARIZATION 

VERTICAL POLARIZATION 

MEAN ERROR 

STD. DEV. 

MEAN ERROR 

STD. DEV. 

13 

-3.87 

11.21 

-4.19 

9.42 

14 

1.49 

12.78 

2.87 

5.45 

16 

3.19 

10.77 

3.81 

6.85 

17 

-2.21 

8.45 

-3.34 

4.57 

18 

-0.81 

8.90 

-1.23 

5.42 

19 

-1.55 

11.95 

-0.75 

7.43 




























(3) Estimate direction using coefficients from (2) and a 0 measured 
with antenna 1. 

(4) Estimate direction using coefficients from (2) and a 0 measured 
with antenna 2. 

(5) Average estimates of direction from (3) and (4) to produce new 
estimate of direction. 

(6) Determine quadrant by comparing the 4 estimates of (5) with data 
from other sources such as: 

(A) concurrent satellite photos 

(B) meteorological analysis and forecasts based on 

observations at other times and places. 

The results show that vertical-polarization is more effective in 
determining both wind speed and direction than horizontal-polarization. 


3.0 COMPARISON OK SOME SCATTERING THEORIES WITH RECENT SCATTERQMETER 
MEASUREMENTS 


A. Introduction 

In recent years there has been general agreement that the scattering 
properties of the sea surface can be explained reasonably well In terms 
of a two-scale roughness model. It is understood in the two-scale model 
that the small waves are assumed to satisfy the small perturbation assump- 
tions and that the large scale waves are to satisfy the Krichhoff approxi- 
mation. This means that only certain high and low frequency portions of 
the entire sea spectrum are included as significant contributors to sea 
scatter. The basic approach is to use the first order result of the small 
perturbation method to compute the scattering coefficient due to small 
scale waves and to account for the tilting effect of the large scale 
waves by averaging this scattering coefficient over the slope distribution 
of the large scale waves. The tfTtfng process not only produces an 
averaging effect but also causes depolarization of the incident 
electromagnetic field in the local reference frame leading to additional 
higher order terms (Fung, J 976) . However, these higher order terms are 
unimportant when dealing with small sea slopes. As a result most authors 
ignore these higher order terms. It is known that the first order per- 
turbation approach invariably produces a larger scattering coefficient 
for the vertical polarization than that for the horizontal polarization 
near grazing incidence. This fact remains unchanged after averaging. 

Since some recent near-grazing observations (Long, 1 974) indicate that 
there are cases where the converse is true, special care must be exercised 
when applying such a theory to explain data at incident angles greater than 
80°. As Wright (1968) pointed out, this type of theory is valid only in 
the mid-angular range approximately 30° to 80°. 

Sea scatter theories using the two-scale model (Wright, 1968, Chan 
and Fung, 1977, Bass, et al., 1968, Wu and Fung, 1972, Long, 197^, Valen- 
zuela, 1967) may be further categorized into two types: (a) those which 

use Phillips* (I966) sea spectrum together with a slope distribution of the 
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large scale waves defined with respect to the horizontal plane, and (b) 
those which use a growing sea spectrum (Pierson, 1 975 > Mutuyasu and Honda, 

1 974, Sutherland, 1567) and a slope distribution of the large scale waves 
defined with respect to the plane normal to the look direction, in the 
next section we shall briefly describe these two types of theory follow- 
ing Wright for type (a) and Chan and Fung for type (b). The predictions of 
these theories are compared with experimental observations in Section C and 
conclusions are given in Section D. 

B. The Two-Types of Two-Scale Models 

Wrfght*s Sea Model . Many authors (Wright, 1968, Chan and Fung, 1977 
and Valenzuela, 1967) have shown that for a slightly rough surface which 
satisfies the small perturbation assumption, the backscatteri ng coefficient 
is given by 

a (0,<i>) = Sk^a^a | 2 W(e ,$) (3-1) 

P P PH 

where for horizontal polarization, p = h and 
°hh ■ cos2eR h 

R. is the Fresnel reflection coefficient for horizontal polarization defined 
n 

in terms of the electric field. FoP vertical polarization, p = v and 

a w = R v cos 2 0 + (k* 2 - k 2 )T v 2 sin Z e/(2k' 2 ) 

where R y and T w are the Fresnel reflection and transmission coefficients 
defined in terms of the magnetic field for vertical polarization. In a , 

r r 

k is the wave number in air; k' is the wave number in sea water and 0 is 
the angle of incidence. W (8 ) Is the normalized anisotropic sea spectrum 
and a is the variance of the small scale sea waves. To include the tilt- 
ing effect of the large scale waves, let j should first be modified by 

rr 

including higher order correction terms due to tilting and then (3-1) is 
averaged with respect to the slope distribution functions of the large 
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scale waves recognizing that 0 in (3-1) is to be replaced by the local 
incidence angle, e 1 , which is a function of the local slopes. Wright 
indicated that satisfactory approximation could be achieved by assuming 
the spectrum to be independent of the tilt angles except as they affect 
the local incidence angle. Letting a be the tilt angle in the plane of 
incidence and 3 be the tilt angle in the orthogonal direction, he 
obtained the backscatter ing coefficient as 


a° pp (e) = 8kV jj |G pp (0,Z x ,Z y )| 2 W(ksin0')P(Z x ,Z y )dZ x d2 y 


(3-2) 


where G vv Sa vv (0,) , Z 7 
G hh 3 + 

Jf.e + z,, s '" 2 |’ W 

P{2 »Z ) o (Sirs 2 ) ^exp 

x y 2 S * 

s 2 - 0.004 + 0.81 x 10"3u 
a 2 W(2ksin9') = 0.585 x 1O“ 2 (2ksin0 ')^ 


and U is the windspeed in m/sec measured 12.5 m above sea level. Note 
that higher order correction terms due to tilting are considered negli- 
gible for vertical polarization so that G yv = a vv ( 0 ’) while for horizontal 
polarization a correction term is added to 0^(0') to form G^. The slope 
distribution, P(Z^,Z^}, is a simplified form reported by Cox and Hunk (1954). 
From the definition of the sea spectrum it is clear that it is independent 

of wind speed. Thus, wind dependence in a° (0) comes only through 

PP 

the slope distribution, p(2 ,Z ). 

Chan and Fung's Model. The model reported by Chan and Fung (1977) 
for the backscattering coefficient is as follows: 


a° (0 ) 
PP 


8k V 



I «pp 1 WfG'^JPgtZx' ,Zy')dZxdZy 
xnt0 


(3-3) 


which is essentially the same as (3-2) except for the definitions of P(Zx* ,Zy') 

and W (e ' ,<{> ) and the fact that higher order corrections to ct are ignored. 

. PP 

tn Chan and Fung (1977) it was shown that the probability of occurrence 

of a given slope varies with the direction of observation. For this reason 
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the* probability distribution function P(Zx',Zy‘) reported by Cox and 
Hunk is projected along the look direction first before it is used in 
the averaging process, i.e. 

P Q (Zx*,Zy') = (1 + Zxtane)P(Zx\Zy') (3-4) 

In (3-4) , P(Xx', Zy') is the slope distribution of the large scale waves 
as viewed at an incidence angle 6 and is defined in the prime coordinates 
(Figure 8) whose x'-axis is parallel to the wind direction. It is assumed 
that the plane of incidence is the XZ plane and that the angle between 
the x-axis and the x'-axis is $ so that an upwind observation occurs 
when $ = 0. In view of Figure 1, the slopes in the primed and the un- 
primed coordinates are related as follows: 

Zx 1 = Zx costf) + Zy sirxfr 
Zy 1 = Zy cos$ - Zx sinji 

The form of the slope distribution of a sea surface, PtZx'jZy*), was 
reported by Cox and Munk (1954) to be 
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c 0,4 ” ®* 2 ^ 






= 0.4 C 22 = 0.12 

= 0.001 - 0.0086U 

= 0.0^ - 0.033U 


- v /(1 + R 0 > 


■ H 

0 U -Q -2 

= {0.003 + 1.92 x 10 5 u)/3.l6 X 10 5 U 


■f. 


0.359 


K 2 s (K)dK 


In the expressions above U is the wind speed in meters per second 
at an altitude of 12.5 m above the sea horizon. The definition of s(K) 
is given in accordance with Pierson (1975) as 


S(K) . 

= s.00. 



K. 

, < 

K < K. 



J 



i- 

-1 
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S,(K) 

3 

“ -qf exp 
K- 3 

f- . 0-7*S 2 , l 
L k 2 u(uJ 4 J 

9 

0 • 

< K « 

:K,=- 

u * 2 

s 2 (k) 

= =K,- ,/2 

k-5/2. 


K 1 

< K 

<k 2 = 

0.359 

s 3 (k) 

‘ V9 

(K/K 3 ) q , 


K 2 

< K 

<K r 

0.942 

s 4 (k) 

= 0.875(2ir) Pj_1 

3gK 2 /l3. 1769 
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(gK + gK 3 /13.1769) tp 1 * 1}/2 
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s 5 (k) 

= 1 .473 x 
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K 4 
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can be found numerically by setting S^K^) equal to S^(K^). 

IL “ friction velocity, Uj. > U,~T 
. f0 - ~rt 

K m = (13.1769) 

q = log 1 q ES 2 (K 2 )/S I( (K 3 ) 3/log, 0 (K 2 /i< 3 ) 

P] = 5-0 - log, 0 U A 

z = 0 . 684/lL + 4.28 x 10‘ 5 uJ - 4 .I 3 x lo“ 2 

Q ** <» 

U(U.J = (U.,/0.4) in (z/z ) cm/sec 

^ o 

a = 05 x 10 _J 

2 

g ~ 980 cm/sec 
lL m = 12 Cifl / sec 
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The Important characteristics to be noted In this spectral model are: 
(0 the spectrum grows with the wind; (ii) in thui capillary region 
defined approximately by S^CK) the larger the K number the faster is 
the growth (i.e. faster growth of the sea spectrum occurs at larger 
incident angles and higher frequencies); and (iii) this model is 
valid to about 38 knots (Pierson and Stacy, 1973). In accordance with 
Mitsuyasu and Honda, their model can be valid to velocity as high as 
33m/sec at 10 meters above the sea horizon. Thus, the theory may be 
valid to higher wind speed than the 38 knots when operating frequency 
is such that only their portion of the sea spectrum is the significant 
contributor. 

The relationship between S(K) and W(0 , ,4i) can be shown to be 
(Chan and Fung, 1977) 


ct Z W(K,$) - f£2- [1 + 2 cos2i]>3 


(3-5) 


where K = 2kstn ' and R is the ratio of the crosswind slope variance 
to the upwind slope variance. Note that R Q is the same as R except 
that it is for the large scale waves only. It Is believed that R 
should be somewhat less than R 0 . However, since the exact expression 
is not known, the theory is illustrated using Cox and Hunk's clean 
sea model for R. Upon substituting (3~*0 , (3“5) into (3 _ 3) we obtain 


V e -*> - 



-k 3 |„ | : 

IT 1 PP 1 


[1 + 


2(i-R) 

1+R 


cos2 ] 


S(2ksin6') 

sinB 1 


-eotO 


(1 + Z^tanQ) P (Zx* , Zy') dZ x dZ y 


(3-6) 
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Equation (3"2) and (3-6) represent the two different types of sea scatter 
theory considered In this paper. Note that the lower limit for Zx in 
(3-6) Is -cot8 instead of minus infinity. This is necessary in order 
to avoid looking at the back side of a facet. As a result, the 
probability density function should be renormalized with respect to 
this range of limit which varies with the incident angle. 

C. Comparison with Sea Scatter Measurements 


1 . Comparison of g°(8,<$) versus Azimuth 

The a°(0,$) versus <{> data reported by Jones and Schroeder (1977) 

are all at 30° incidence angle. Comparisons between Chan and Fung's 

theory and these data are shown in Figure 9a and 9b for the vertically 

and horizontally polarized cases respectively. Wright’s theory is not 

shown since it does not include azimuth dependence, in making these 

comparisons, the absolute levels of a 0 and a?, at different wind 
• vv hh 

speeds are disregarded, it is seen that satisfactory agreement is 
obtained for the variations of o° versus $ for wind speeds at 6.5 
m/sec and 15 m/sec and for both polarizations. As will be shown in 
0°(0,$} versus wind speed comparisons in Figure 10, it is rather mean- 
ingless in most cases to compare the absolute levels of o°(0,<ji) at 
a single incident angle, polarization, and wind speed. This comes 
about because of the unavoidable scatter which is always associated 
with limited amount of data. 

2. Comparison of g°(8,ji) Versus Wind Speed 

Figure 10a, b, and c, show the comparisons between ff° data re- 
ported by Jones and Schroeder (1977) and 4^ vv (0,<J>) computed using 
(3 _ 2) and (3-6) versus wind speed at 19.5 m above the sea horizon 
for downwind and crosswind cases respectively at incidence angles 
30°, *jQ°, and 50°. In these comparisons the computed absolute levels 
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Figure 9a: Theoretical and observed values of o° plotted versus the azimuth angle. 
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of vv (0 are retained. The range of wind speeds considered is ap- 
proximately between 6.2 m/sec and 15 m/sec. It is noted that the reason 
why is much higher in level than crj^(e,$) is the inherent char- 

acteristics of the first order perturbation theory. While the MFE data 
are limited, they seem to agree with the predicted theoretical level se- 
paration. Also, uncorrected data. Figures l la and b, taken from Skylab 
5L-2 and SL-3 at ^3° and 50° do show a clear level separation of the pre- 
dicted order of magnitude. Note that data selected from Skylab are those 
taken during Hurricane Ava and Tropical Storm Christine, where a continuous 
measurement of o° versus wind speed, was available. Despite clear differences 
at particular incident angles and wind speeds, a general agreement in absolute 
levels in wind speed trends and in incident angle variations is observed. 

Also included in Figures !0a and 10b are the theoretical predictions 
from Wright’s theory. It is seen that g° exhibits much less dependence 
on wind speed. This appears to be due primarily to the use of Phillips’ 
sea spectrum model which is independent of the wind speed. Recent studies 
on the hydrodynamic effects of the long and the short wave interactions 
may serve as another explanation for the observed increase in o°. How- 
ever, since non-linear Interaction is of higher order, It does not seem 
likely that a full explanation can be made in terms of this effect alone. 

In view of the fact that the measured a° values do increase significantly 
with the wind, it appears that the two-scale scatter theory cannot provide 
correct predictions unless the sea spectrum is allowed to increase with 
the wind speed, 

D. Conclusions 


When an appropriate sea spectrum and sea-surface slope distribution 
for the large scale waves are Incorporated, the two-scale scatter model 
appears to provide satisfactory explanations of the variation of the 
scattering coefficients versus polarization, azimuth, and wind speed. 
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k.O CONCLUSIONS AND RECOMMENDATIONS 


The 1975 JONSV/AP circle flight data appears to be adequately modelled 
by 


a°(U,0 ,0 W >= A o fU,e) + Aj (U,8)cos((ji w ) + A 2 (U,e)cos (2^ w ) (2-3a) 


The coefficients of this model, determined through regression analysis 
were shown to be related to wind speed by 


A n (u,e) 


Y n f0) 

p n (e) u n 


(2-k) 


However, previous circle flight data sets exhibit an assymetry not 
accomodated by this model, Afarani (1975). Explanations for these 
observations Include both experimental problems and the possibility that 
they are caused by the effect of underlying wave or swell trains traveling 
in a direction different from the wind direction. The difficulty of 
making such measurements from an aircraft cannot be overestimated. Not 
only must the aircraft maintain a constant altitude during the circle, but 
either the wind conditions must be the same at all points in the five- 
mile diameter circle or some kind of correction must be made for differences 
that can only be inferred because wind measurements are made^at a point 
near but not in the circle. Consequently, measurements from a tower 
seem in order since the wind can be measured on the tower and no obser- 
vation point is far enough away to be in a different wind field. 

The feasibility of the WINDVEC inversion algorithm to determine the 
wind speed and direction from the orthogonal scattering measurements was 
demonstrated. In this algorithm, an estimate of the wind speed is 
obtained from the average of the radar cross section by the power-law 
relationship. 


U = a(0) 



(2-7) 


Using this estimated wind speed values for the A coefficients in (2-3a) 

n 

can be obtaind. Clearly, the only remaining unknown in this expression 
Is the wind direction, which can be solved for by 


* 


w 


cos 


-I 
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(2-9) 


This yields two estimates for the wind direction, one tn the first 
quadrant and the other in the second. These values are mirrored to obtain 
values in the third and fourth quadrant, respectively. The estimates of 
wind direction from the fore and aft beam are corrected to North and then 
averaged to improve the wind direction estimate. Vertical polarization 
estimates of both wind speed and direction are considerably better than those 
for horizontal polarization. 

Satisfactory application of this inversion to SASS (SeaSat Active 
Scatterometer System) will require a larger data bank with a wider range 
of incidence angles. At present sufficient good data exists only for 0 = 40°. 
Consequently, before SeaSat launch there is a need for more experimentation, 
either tower or aircraft, so that optimum advantage can be obtained from its data. 

Also, comparison between different scattering theories and scatterometer 
measurements have shown different results. First, better agreement is 
obtained when the slope distribution reported by Cox and Hunk is projected 
along the look direction before it is used in the averaging process. 

Second, use of a sea spectrum which grows with wind speed also yields more 
satisfactory results. 
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Wind Direction Error: 

(PH I ERR) vs 4 

A. Direction errors using antenna 1. 

B. Direction error using antenna 2. 

Wind Direction Error: 
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